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I  Introduction 


Simulation  models  of  two  representative  STOL  aircraft 
have  been  generated.  The  models  are  documented  in  this 
report. 

The  computer  simulation  is  to  be  used  as  a  tool  in  the 
development  of  STOL  terminal  area  guidance  and  navigation 
systems . 

This  intended  use  has  determined  the  form  of  the  simulation: 
The  aircraft  are  described  by  means  of  non-linear  equations  that 
will  accomodate  gross  changes  in  angle  of  attack,  pitch  angle,, 
flight  path  angle,  velocity,  and  power  setting.  Aircraft  motions 
ii:  response  to  control  inputs  and  external  disturbances  are 
related  to  Earth-fixed  coordinates.  The  aquations  are  programmed 
to  run  in  "real  time”  so  that  they  can  be  used  in  conjunction 
with  a  manned  cockpit  simulator.  Provisions  are  made  for  pilot 
control  inputs  to  the  simulation,  and  conventional  panel  display 
parameters  are  generated. 

The  aircraft  which  are  modeled  -  the  DHC  ”Twin  Otter”  and 
the  DHC  "Buffalo"  -  are  described  in  Figures  i  and  2,  respectively. 
They  were  selected  as  representative  light  and  medium  propeller- 
driven  STOL  transports.  Their  selection  does  not  imply  that  there 
are  not  other  STOL  aircraft  representative  of  these  classes . 
Similarly,  the  material  contained  in  this  report  should  not  be 
used  as  the  basis  for  an  evaluation  of  the  flying  qualities  of 
the  "Buffalo”  or  "Twin  Otter"  or  of  the  suitability  of  these 
aircraft  for  any  specific  mission. 


The  aircraft  are  modeled  only  to  the  extent  necessary 
to  yield  a  representative  vehicle  model  controllable  by  a 
guidance  or  navigation  system.  Certain  simplifying  assumptions  - 
specified  in  the  following  sections  -  are  made.  These 
assumptions  are  justified  for  the  present  model  application 
but  may  render  the  model  unsuitable  for  other  possible  applica¬ 
tions  . 

The  simulation  is  described  in  detail  in  the  following 
sections  of  this  report.  In  Section  II,  all  required  equations 
are  developed.  Section  III  tabulates  numerical  values  to  be 
used  in  these  equations  for  the  "Buffalo”  and  "Twin  Otter" . 

The  simulation  program  is  presented  in  Section  IV.  A 
listing  of  all  computer  statements  is  included.  Finally,  in 
Section  V,  representative  simulation  results  are  shown.  These 
results  demonstrate  that  the  simulation  is  an  adequate  representa¬ 
tion  of  the  two  STOL  aircraft. 
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II  Description  of  Mathematical  Model 


The  mathematical  model  consists  of  all  equations  required 
to  describe  the  motions  of  the  aircraft  in  space  resulting 
from  external  disturbances,  control  inputs,  and  the  aircraft's 
aerodynamic  characteristics.  These  equations  are  presented 
in  this  Section.  First,  however,  it  is  necessary  to  define 
the  reference  coordinate  frames  to  be  used. 

IIA  Definition  of  Reference  Coordinate  Frames 

Reference  coordinate  frames  to  be  used  in  this  analysis  are 
defined  in  this  section.  Insofar  as  possible,  axis  systems  have 
been  defined  so  that  senses  of  rotation  and  translation  are 
similar  for  small  rotations.  Positive  force,  moment,  and  motion 
vector  components  are  defined  to  be  in  the  positive  sense  of 
the  axis.  To  the  largest  extent  possible,  the  symbols  and 
conventions  used  are  consistent  with  those  in  common  usage  in 
the  guidance  and  control  fields  and  with  those  used  by  NASA  for 
aircraft  stability  and  control  work. 

The  Earth  Local-Vertical  Frame  (L)  is  a  local  geographic 
frame.  Its  origin  is  fixed  at  a  point  on  the  Earth's  surface  with 
Z,  along  the  vertical  defined  by  the  local  gravity  vector 

ij 

(positive  downward)  ,  Xj^  parallel  to  geographic  North  (positive 
to  the  North)  ,  and  parallel  to  geographic  East  (positive 
to  the  East) . 

The  Aircraft  Body  Coordinate  Frame  (A)  is  fixed  to  the 
aircraft  and  r'*-^  as  and  translates  with  the  aircraft.  Its 
origin  is  the  center  of  mass  of  the  aircraft.  The  axis 
is  chosen  in  a  forward  direction  in  the  plane  of  symmetry  that 
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i 

1 

1 

i 
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is  parallel  to  the  initial  or  equilibrium  direction  of  the 
remote  wind.  Thus  the  A-frame  axes,  by  the  commonly  accepted 
definition,  are  "stability  axes".  Because  the  axis  is 
initially  aligned  with  the  remote  wind,  the  initial  angle  of 
attack  a(o)  =  is  zero.  (In  this  report,  0  and  a  when  not 
subscripted  to  indicate  reference  frame ,  are  assumed  to  be 
referenced  to  the  A-frame.  Further,  since  in  the  simulation 
documented  in  this  report  the  aircraft  is  placed  in  equilibrium 
at  t  =  0,  "equilibrium"  and  "initia"' "  conditions  are  equivalent.) 
The  Y  axis  is  normal  to  the  aircraft's  plane  of  symmetry 
(positive  to  the  right) ,  and  the  Z  axis  is  in  the  plane  of 
symmetry  (positive  downward)  and  orthogonal  to  the  and  Y^ 
axes.  The  A-frame  is  related  to  the  L-frame  (and  to  the  next- 
defined  C-frame)  in  Figure  3. 

The  Earth-Aircraft  Control  Coordinate  Frame  (C)  is  also 
centered  at  the  center  of  mass  of  the  aircraft.  The  axis 
is  aligned  with  the  local  gravity  vector  (positive  downward) 
and  is  therefore  parallel  to  the  Zj^  axis.  The  X^  axis  is 
the  intersection  of  the  horizontal  plane  with  the  vertical 
plane  containing  the  X  axis .  The  Y  axis  completes  the 
orthogonal  right-hand  system.  The  C-frame  is  an  intermediate 
frame  needed  to  define  the  Euler  angles  describing  the  relation¬ 
ship  between  the  Earth  local-vertical  (L)  frame  and  the  Aircraft 
body  (A)  frame.  In  their  order  of  rotation  (which  must  be 
preserved)  the  Euler  angles  are  defined  as: 
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1.  Heading  (W) :  angle  of  rotation  about  Z 

‘  from  X  to  X^/  ■  I  I  i 

L 

2.  Pitch  (0)':  angle  of,  rotation  about  Y_ 
from  X  to  X, ;  ' 

3.  Roll  ($)  :  an'gle  of  rotation  about  X 

from  to  Y,  ^ 

C  A  :•  I 

•  * 

*  I  : 

These  Euler  angle  rotations  are  shown  in  Figure  3. 


The  Body  Reference  Coordinate  Frame  (B)  js. introduced  and 

I 

defined  in  this  report  primarily  'to  clarify  the  definition  of 

.  *  ' 

trim  angle  of  attack.  Like  the  A-frame,  this  frame  is  fixed  , 

! 

to,  and  translates  and  rotate's  with,  the  aircraft  and  has,  as 

I 

I  I 

its  origin  the  center  of  mass  of  the  aircraft.  The  X^  axis, 

however,  is  fixed ’in  a  forward  direction  in  the  plhne  of 

symmetry  parallel  to  a  fuselage  waterlinfe  or  datum  line..  The 

'  1 

X_  axis  IS  displaced  from, the  X.  axis  by  the  angle  a'  .  The 
“  ■  ,  •  I  A  I  •  •  "O  I 

Yg  axis  coincides  with  the  Y^  ax’is,  and  the  iZg  axis  (positive  ' 
downward)  forms  an  orthogoi^al  set. 

'  I 

,*  the  trimmed 

angle  of  attack.  ’  It  is.  the  .angle  between  ^:he  initial  (equilibrium) 

*  »  ! 

J  • 

remote  wind  vector  and  the  X_  axis.  Unlike  Oq,  it  has  a  non  -  ’ 

*  I  "  ' 

zero  value.  It  is  evident  from  Figure  4  ‘that 


The  angle  a„  is  sometimes  'cal  led  la^  .  , 
'  trim' 


00  =  0  +  ag 


(D*  ‘ 


^Numbered  equationj  are  mechanized  in  the  simulation.  Other 
equations  are  introduced , as  necessary  for  purposes  of  clarification, 
but  are  not  numbered:. 
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and,  for  equilibrium  levfel  flight,  that 

•  .  *  » 

0T3  “•  CCq  ■  •  '  ‘  ^  I 

'  O  ®o 

I 

^  ■  ' 

,  The  trim  angle  Un  can  be  approximated  in  the  following 

“O  '  ' 

fashion  in 'the  absence  of  wind  tunnel  or  flight.,  test  =data« 

'  '  I 

'  Assumihg  a*  constant  j»'  ...raft  lift  durve  slope,  a,  :sketch 
.  '  • 
'  ‘  I 

the  aircraft's  lift  curve: 


a. 


B 


or,  at  equiljibrium. 


C,  =a  (a„  ) 


Next,  assume  that  wing  incidence  has  been  chosen  by  the  aircraft 
manufacturer  to  produce  a  .level  fuselage  attitude  '{tto  =  0)  when, 

I  1  .  "o  • 

'the  aircraft' is  in  flight  at  "Economy  Cruise  Speed"  at  10000  ft 

•  I  • 

I  ^ 

and  at  an  arbitrarily  -  'chosen  average  gross  weight.  Usings  the 

I  ' 

relation  W  _  =  Cr  q  S,  calculate  the  .lift,  coefficient  at  the 
cr  .  L,  II. 


cr 


•cr 


flight  condition.  The  ?  ^jie  of  'attack  for  zero  lift  c^n  then 


i  , 


be  calculated  from  the  .»bove  equation  as 


Cl, 


a 


cr 


B 

OL 
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The  same  equation  can  be  manipulated  to  give  an  expression 


for  the  trim  angle  Ug  at  any  other  trim  lift  coefficient: 

Cl 

^o 

“b.  =  —  +  “B 

w  f\r 


(In  Appendix  B  of  Reference  1#  C_  was  estimated  to  be  .44 

^cr 

for  the  "Buffalo”  and  .48  for  the  "Twin  Otter".  For  both 
aircraft/  a  =  5.2/rad/  so 


cto  =  -.035  =  “4.8®  (Buffalo) 

OL 

=  -.092  =  -5.3®  (Twin  Otter) 
These  values  are  used  in  this  report.) 


IIB  Velocity  Resolutions 

Use  must  be  made  of  the  above-defined  Euler  angles  to 
relate  a  vector  quantity  in  the  A-frame  to  its  components  in 
the  L-frame  and  vice  versa.  In  general/  a  vector  R  can  be 


resolved  into  its  A-frame  or  L-frame  components: 


R  =  By  ^  + 


.  j  .  +  R  k 

a^a  ^ 


°  \  ^  "zl  ’'l 

Where  i /  j ,  and  k  are  unit  vectors  in  the  indicated  frames . 

L-frame  components  of  R  can  be  expressed  in  terms  of 
A-frame  components  of  R  and  the  Euler  angles: 


“  7  - 


where  =  cos  T  cos  0  (3) 

®12  ~  ®  ^  ”* 

=  cos  H*  sin  0  cos  t  +  sin  '{'  sin  ^  (5) 

02^  =  sin  'H  cos  0  (6) 

022  =  sin  "J*  sin  0  sin  $  +  cos  ’?  cos  *  (7) 

022  =  sin  '1'  sin  0  cos  t  -  cos  ’{'  sin  'I'  (8) 

B^j^  =  -sin  0  (9) 

0^2  =  cos  0  sin  4>  (10) 

0^^  =  cos  0  cos  <1>  (11) 
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Conversely,  A-frame  components  of  any  vector  ^  can  be  expressed 
in  terms  of  L-frame  components: 


-8- 


Thus,  in  the  simulation,  the  A-frame  components  of  aircraft 

velocity  (x  =  u,  :^  =  V,  and  4.  =  W)  are  computed  and  used 

A  A  A 

to  obtain  velocity  components  with  respect  to  the  ground;  , 


X  =  B 
L 


Yr.  =  B 


u  + 

B 

V  + 

W 

(fps) 

(12) 

11 

12 

13 

'21 

®22 

V  + 

®23 

W 

(fps) 

(13) 

■®31  “  • 

“®32 

V  - 

®33 

w 

(fps) 

(14) 

I 


r 


I 


lie  Provisions  for  Atmospheric  Disturbemces  (Winds) 

Winds  are  input  into  the  simulation  in  the  L-frame. 

Components  are  (positive  North)  ,  y,^  (positive  East)  , 

L  L 

and  zw  (positive  downward) .  The  winds  are  resolved  into 
L 

A-frame  components  in  equations  15-17  in  order  to  compute 
airspeed  components ; 


°  -  [Bn  \  y„^  + 

• 

V=V-fB  *  +B^V  +B  z] 

w  ''  22  32 

W  +  B  y  +  B  z  ] 

w  13  W^  23  \ 


(fps) 

(fps) 

(fps) 


(15) 

(16) 
(17) 


Material  contained  in  this  report  is  sufficient  to  allow 

introduction  of  steady  state  wind  components.  The  desired 

•  •  • 

winds  are  simply  input  as  x  ,  y„  ,  and  z  .  The  report  does 

L 

not  document  wind  gust  or  wind  shear  models.  However,  these 
models,  when  developed,  can  be  readily  incorporated  into  the 
simulation  with  only  minor  modifications  to  the  program  being 
required. 
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IIP  Airframe  Equations  of  Motion 


In  Reference  1,  general  6  degree  of  freedom  airframe 
equations  of  motion  were  developed  as 

m  [U+QW-RV+gsin0]  =X  (longitudinal  force) 

m  [V  +  RU  -  PW  -  g  cos  0  sin  $]  =  Y  (side  force) 

m  [W  +  PV  -  QU  -  g  cos  0  cos  $]  =  Z  (normal  force) 


P  +  (I  -  I  )  QR 
X  z  y 

I  Q  +  (I,,  -  I,)  RP 

y  X  z 

I  R  +  (I  -  I  )  PQ 
z  y  X 


-  J 

(R  +  PQ)  =  L 

xz 

- 

S 

II 

1 

CM 

xz 

-  J 

(P  -  QR)  =  N 

xz 

(rolling  moment) 
(pitching  moment) 
(yawing  moment) 


where  the  body-axis  angular  rates  P,  Q,  and  R,  can  be  used  to 
obtain  Euler  angle  rates  according  to  the  equations 


W  3  Q  ^  ^  T>  cos  4* 

cos  0  COS  0 

0  =  Q  COS  4>  -  R  sin  * 
$  =  P  +  y  sin  0 


(rad/sec) 

(18) 

(rad/sec) 

(19) 

(rad/sec) 

(20) 

These  nine  equations,  together  with  equations  12-14,  provide  an 
almost  exact  description  of  the  motions  of  an  aircraft  operating 
near  the  Earth's  surface.  They  involve,  as  shown  in  Reference  1, 
only  four  assumptions : 

1 .  Aircraft  mass  is  constant 

2 .  The  Earth  can  be  considered  an  inertial  frame 

3.  The  aircraft  is  a  rigid  body 

4 .  The  aircraft  is  symmetrical  about  its 
X  -  z  plane. 
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For  purposes  of  this  simulation,  the  above  6  rigid  body 


airframe  equations  have  been  approximated  as 

.  2 

U  =  RV  -  QW  -  g  sin  0  +  X/m  (ft/sec  )  (21) 

V  =  PW  -  RU  +  g  cos  0  sin  4  +  Y/m  (ft/sec^)  (22) 

•  2 

W  =  QU  -  PV  +  g  cos  0  cos  4  +  Z/m  (ft/sec  )  (23) 

P  =  L/I^  (rad/sec^)  (24) 

Q  =  M/I  (rad/sec^)  (25) 

.  ^  2 
R  =  N/I^  (rad/sec  )  (26) 

z 


The  omitted  terms  in  the  moment  equations  involve  either  products 

of  angular  velocities  (e.g.  QR)  felt  to  be  small  compared  with 

other  equation  terms,  or  terms  containing  J  which  will  be 

xz 

neglected.  Experience  has  shown  that,  for  purposes  of  this 
simulation,  these  terms  can  be  omitted  with  negligible  effect 
on  results. 

The  terms  X,  Y,  Z,  L,  M,  and  N  of  equations  2i  >-  26  represent 
the  aerodynamic  forces  and  moments  acting  on  the  aircraft.  The 
lateral  terms  (Y,  L,  N)  will  be  expressed  in  a  quasi-linear  form 
(as  in  Reference  1),  but  the  longitudinal  forces  and  moment  (X,Z,M) 
must  be  non-linear  in  order  to  permit  large  excursions  in  forward 
velocity . 

The  longitudinal  aerodynamic  force  terms  are,  from  the  sketch, 

X  =  T  -  D  cos  a  +  L  sin  a  (lbs)  (27) 

Z  -  -  (Laos  a  +  D  sin  a)  (lbs)  (28) 
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The  terms  X  ,  Z_,  Z»,  and  Z^  have  been  neglected  in  this 
q  q  w'  Og 

analysis  because  of  their  small  contribution  to  the  overall 
forces . 

It  is  also  assumed  that  all  thrust  forces  act  along  the  X 

A 

axis.  Thus  moment  effects  of  thrust  changes  are  neglected,  as  are 
forces  and  moments  produced  by  special  lift  devices  operating 
within  or  outside  of  the  propeller  slipstream.  These  effects 
are  neglected  because  the  airframe  data  required  to  model  them 
are  not  available. 

Equations  27  and  28  are  solved  .(as  are  the  other  simulation 
equations)  once  every  computer  iteration  cycle.  Thrust,  drag, 
and  lift  force  components  are  summed  to  produce  resultant  X  and 
Z  forces  acting  on  the  aircraft. 

Expressions  for  the  total  thrust,  lift,  and  drag  forces 
are  next  developed. 

Thrust  is  computed  from  an  empirically-derived  expression 
(developed  in  the  appendix)  which  accounts  for  the  effects  of 
altitude  h,  airspeed  V  ,  and  throttle  setting 


Nsaamoi 


i 


\ 


a  T 


T  = 


static 


1  \  ^  ?T2 


T  •  5 


■(lbs) 


■  (29)  • 


where  0  5  5  1  1.0, 


tl 

-h/h  ^ 

a  =  e  atm. 

1 

(-) 

.(30) 

1 

4 

V„  =  [U, +  V  ^  +  W  2] 
R  w  w  w 

1/2 

1 

1  1 

(fps) 

1 

! 

(3lj 

1  Lift  and  drag 

are  calculated  from  the 

'standard 

relationships : 

L  =  Cl  qs 

(lbs')  1 

(32) 

D  =  Cjj  qS 

■  (lbs)  ' 

(33) 

1  where 

1 

!» 

i 

Cl  =  Cl^  +  aa 

» 

1 

<-)  ■  , 

1 

(34) 

1 

2  , 

C^  =  Cn  +  Cr  /freAR 

D  Df  L 

('") 

f 

(35) 

J 

2  R 

•  (Ibs/ft^) 

.  (36)' 

1 

i 

P  =  a  Pq 

i 

(sl/ft^) 

(37). 

1 

1  and 

a  =  tan"^ 

1 

(rad) 

(38) 

The  expression  for  pitching  mom  at-  used  in  the  simulation 


IS 


M  =  qSclCn,^  +  •  (ft/lbs)  (39) 


a 


«  R  -  -  -e 

where  the  coeffi2ients  of  the  variables  are  constants.  The 

term  is  zero  in  this  report,  but  is  included  to  facilitate 
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I 


I 


later  shaping  of  the  trimmed  h-  vs  V_  curve.  To  do  this,  C_ 

■  ,  e  K  III. 

I  '  ^ 

would  be  lAade  a  function  of  V  . 

R 

Rate  of  change  with  time  of  angle  of  attack  is  obtained  by 
.  <  ' 

jdiffere’ntiating  equation  38;  ’ 

‘  ,  Wi 

-  at  0  >  ■ 

'  .  '  '  . 

:  U„  W  -  W  U 

i  w  W  w  w 


u, 


v; 


2  +  W  2 
.  w 


If  the  approximation  is.  made  that  U  =  U  ‘and  W  =:  W  ,  the  above 
;  I  w  w 

expression  can  be  manipulated  to  produce 


'  W  I  ■  ’2  , 

.  I.  w  .  “ 

a  =  '(W  -  U) 


U. 


U 


(rad/sec) 


w 


w 


(40) 


I  1  which  is  the  expression  used  in  the  'simulation. 

'  The  lateral  force  (Y)  and  moments  (L  and  N)  are 

I  (  !  ^  * 

's  developed  in  conventional  linearized  form  (as  in  Reference  1) 

except  that  total  variables  are  used  rather  than  perturbation 

I  •  ■  . 

i  values,  and  that  coefficients  of  the  lateral  variables  are  made 

i  I 

if unctions  of  lift 'and  drag  coefficient,  airspeed,  and  dynamic 

I  I  • 

pressure,  all  of  which  are  determined  by  solution  of  the 
longitudinal  equations . 

:The  lateral  force  and  moment  expressions  used  in  the 
simulation  are;  '  ,  i  i 


Y  = 

Y  V  +  Y  R,+  Y  P 

V !  w  r  p 

1 

(lbs) 

(41) 

l'  = 

\  Yw  ^  ^ 

I'Sa'  'a 

(ft-lbs) 

(42) 

N  = 

Vw  ^  Np  P  + 

(ft-lbs) 

(43) 
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The  terms  Yg  ,  Lg  ,  and  Ng  ,  sometimes  included  in  the  lateral 
r  r  a 

equations/  have  been  omitted  in  the  present  analysis  because  of 
their  negligible  effects , 

The  coefficients  of  these  equations  are 


=  -  D  V  SC 

2  P  ® 

(Ibs/fps) 

(44) 

=  ^  p  Sb  C 

(Ibs/lU) 

(45) 

4  R  Yj. 

sec 

=  7  p  V„  Sb 

(Ibs/EM, 

(4e) 

4  R  yp 

sec 

=  1  p  Vjj  Sb  c. 

(ft-lbs/fps) 

(47) 

=  i  p  V  Sb^  0, 

4  R 

(ft-lbs/||f) 

(48) 

=  C,  +  0^4 

"'fin 

(-) 

(49) 

1  2 

,rad. 

(50) 

,=IPV^Sb 

i  =  q  sb  Cji 

(ft-lbs/rad) 

(51) 

a  ®a 

,  =  I  P  \  Sb  c„ 

(ft-lbs/fps) 

(52) 

6 

1  2 
=  7  p  V  Sb  C„ 

(ft-lbs/^) 

(53) 

•  4  R 

sec 

‘r 

r  wing 

(-) 

rad 

(ft-lb/— ) 
sec 

(54) 

1  2 
=  7  p  V^Sb'^  C„ 

,  4  P  R  n 

(55) 

_a_ 

1  “  ^n  "4  ^  ”  ttAR 

P  %IN 

(-) 

(56) 

5  =  q  Sb 

(ft-lbs/rad) 

(57) 
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The  equation  for  sideslip  angle  is 


3  =  tan"^  (rad)  (58) 

Uw 


Linear  and  angular  rates  are  integrated  to  produce  the  required 
linear  and  angular  displacements.  Initial  values  of  displacements 
are  provided  for  where  necessary: 


u  = 

t  • 

u(o)  +  /  u 

dt 

(fps) 

(59) 

V  = 

o 

.t  ' 

v(0)  +  1  V 

dt 

(fps) 

(60) 

w  = 

o 

t  • 

w(0)  +  /  w 

dt 

(fps) 

(6)) 

p  = 

o 

• 

J  P  dt 

(rad/sec) 

(62) 

Q  = 

0 

t  • 

/  Q  dt 

(rad/sec) 

(63) 

R  = 

o 

R  dt 

(rad/sec) 

(64) 

y  = 

o 

-t  . 

/  'F  dt 

(rad) 

(65) 

0  = 

o 

ft  • 

J  0  dt 

(rad) 

(66) 

$  = 

o 

ft  * 

/  $  dt 

(rad) 

(67) 

X 

o 

ft  • 

=  J  X  dt 

0  ^ 
ft  . 

=  J  ■  y  dt 

0  ^ 

-  2_  =  h(0) 

Xi 

(ft) 

(68) 

L 

h  = 

+  /  h  dt 
o 

(69) 

(70) 
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HE  Definition  of  Required  Display  Quantities 

Provisions  are  made  in  the  simulation  for  displaying 
parameters  that  are  commonly  available  on  a  cockpit  instrument 
panel.  These  parameters  are  tabulated  here  (and  are  defined  if 
they  have  not  been  previously  defined) : 


Indicated  Airspeed  IAS  =  . 

(mph) 

Altimeter  Output  h 

(ft) 

Directional  Gyro  Output  57.3 

(deg) 

Pitch  Attitude  Gyro  Output  57.3  0^ 

(deg) 

Roll  Attitude  Gyro  Output  57.3  $ 

(deg) 

Rate  of  Climb  Indicator  Output  li/60 

(fpm) 

Turn  Rate  Indicator  Output  57.3  R 

(deg/sec) 

Slip  Indicator  Output 

.g  cos  0  sin  $  -  V  -  RU  +  PW. 
g  cos  0  cos  $  -  ^  -  PV  +  QU 

(rad) 
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Ill  Tabulation  of  Numerical  Data  for  “Buffalo”  and  "Twin  Otter” 
Numerical  data  for  the  two  aircraft  to  be  modeled  are  tabu- 


m. 


I 


lated  in  this  section.  Unless  otherwise  indicated,  the  values 
have  been  taken  from  Reference  1.  It  should  be  recognized  that 
stability  derivative  values  tabulated  here  are  not  based  on  wind 
tunnel  or  flight  test  results,  but  have  been  generated  using 
analytical  expressions  presented  in  Reference  1. 


Parameter  Value 


Buffalo 

Twin  Otter 

a,rad“^ 

5.2 

5.2 

AR 

9.75 

10 

b,ft 

96 

65 

c,ft 

10.1 

6.5 

.032 

.039 

Ap 

Cd 

.030 

.035 

^mt 

0 

0 

-35.6 

-24.6 

mq 

C 

-.78 

-.78 

““a 

‘“m. 

-6.05 

-6.15 

a 

C 

2.12 

1.73 

^6e 

Cp. 

-.53 

-.53 
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Parameter  ’  iValue 


, 

'  Buffalo 

1 

Twin  Otter 

-■.125' 

,  \ 

-!i03 

1 

..20  >  . 

1 

.38 

.  .03^ 

.033 

^fin 

1 

"  s 

Cn  ’  ■. 

!p25 

1 

.033 

Pfin 

1  1 

C 

-.169 

-.168 

”r«- 

■^fin 

! 

'  1  >  , 

Cn 

.101 

1 

.121 

i  * 

C„  1 

.107 

.107 

1 

'  ’  1 

r 

1 

,  -.055 

-.085 

p 

• 

)  1 

=yr 

i 

.368  ’ 

.429 

1 

-.362 

-.492 

1 

1 

C^j^,fps“^  (1) 

.00370 

.00.378  : 

1  ^ 

Cip^/fps-^  (1) 

.6.51x10“®  1 

t 

9.07x10"^ 

e 

\ 

.75 

.7,5 

(2) 

32500 

1 

32500 

Ix/Slug-ft2 

273000 

; 

2430'0 

i 

Iy,slug-ft^ 

i 

215000 

220<^0 

1 

l2,slug-ft^ 

1  V 

447000 

,  41000 

Jjjgi  slug-ft^ 

0 

1 

0  ‘ 
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Parameter  '  ' 

f: 

1  Value  '  '  ! 

:  ! 

>  •  Buffalo 

4 

Twin  Otter 

S,ft^ 

'  945 

1 

420  ■ 

'^static 

22400 

5750  ’ 

W , lbs  ■  ' 

1  ' 

40000  , 

.  12000 

Od  ,rad  (3) 

OL  !  = 

1  I 

'  ■  -.085 

-.092  •  ! 

*  i 

PQ,slugs/ft^  ' 

I  . 

.002378 

’  ’  .002378, 

t 

Notes  1.  Pron{  Appendix,  this  report, 

2,  Atmospheric  density  ratio  calculate’d  as 
'  .  0  =  g”h/32500  compares  with  standard 

atmosphere  data  as  foJLlows:  « 


0 

5000 

10000 

15000 

20000. 


1  1 
.8'62 
,738 
.629 

.533 

1 


.858 

.735 

.630 

.540 
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3.  ■  from'  Section  IIA,.  'this  report 


IV  Simulation  Program 

The  equations  of  Section  II  have  been  programmed  for  real¬ 
time  solution  on  an  XDS9300  digital  computer  at  the  TSC  Simulation 
Facility. 

Because  the  simulation  is  a  simple  one,  a  flow  chart  is  not 
presented.  The  program  listing,  together  with  the  discussion 
presented  here,  should  be  sufficient  to  completely  describe  the 
simulation.  The  listing  is  included  in  this  report  as  Table  I. 

IV -A  Interface  with  GAT-1  Cockpit 

Provisions  are  made  to  drive  the  simulation  manually  using 
a  GAT-1  fixed-base  cockpit  modified  for  the  purpose.  Commands 
from  the  cockpit  are: 

Elevator  trim  (ELTRM)  ‘ 

Longitudinal  stick  displacement  (DLE) 

Lateral  stick  displacement  (DLA) 

Rudder  pedal  diaplacement  (DLR) 

Throttle  setting  (THROT) 

The  scaling  voltages  used  are  given  in  Table  I. 

I 

Similarly,  the  display  quantities  presented  at  the  GAT-1 
panel  (listed  in  Section  II-E)  are  scaled  as  shown  in  Table  I. 


IV-B  Definition  of  Initial  Values  of  Variables 

It  is  convenient  to  be  able  to  begin  a  simulation  run  with 
the  aircraft  trimmed  at  a  level  flight  condition.  Accordingly, 
provisions  are  made  in  the  simulation  for  inputting  desired 
initial  conditions,  and  then  for  calculating  required  initial 
values  of  other  parameters  to  produce  a  trimmed  flight  condition. 


Non-zero  initial  values  are  normally  input  for  altitude 
h{0)  and  airspeed  V  (0) .  In  addition,  non-zero  steady  state 
wind  values  can  also  be  specified.  Zero  initial  values  are 
set  in  the  first  computer  iteration  for  these  parameters: 

U,  V,  W,  P,  Q,  R,  ip,  e,  i>,  P,  Q,  R, 
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HTHT  III!  I 


The  last  two  equations  define  required  pilot  inputs  for 
initial  triw.  In  the  simulation,  provision  is  made  for 
inputting  these  trim  values  for  a  specified  length  of  time, 
after  which  the  actual  control  signal  from  the  cockpit  is 
used.  The  magnitude  of  the  delays  are  TMTHR  seconds  for 
throttle  setting  and  TMDLE  seconds  for  elevator  input  6^. 
This  scheme  permits  setting  up  an  inital  trimmed  condition 
without  the  need  for  cockpit  control  manipulation.  it  is 
useful  when,  for  example,  step  response  runs  are  to  be  made. 
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V  Simulation  Results 

Simulation  results  are  presented  in  this  section.  These 
results  are  in  the  form  of  time  responses  to  various  step  control 
inputs . 

The  time  responses  are  presented  in  a  manner  that  permits 
direct  comparison  with  the  linearized  results  generated  in 
Appendix  D  of  Reference  1.  In  general,  agreement  between  the 
two  sets  of  responses  is  very  close. 

It  should  be  noted,  however,  that  Reference  1  and  this  report 
utilize  the  same  analytically-derived  data.  Therefore  agreement 
between  these  two  reports  does  not  in  itself  prove  the  validity 
of  either  set  of  results.  This  proof  can  only  be  obtained  by 
comparing  the  present  results  with  data  obtained  from  some  other 
independent  source.  Unfortunately,  however,  specific  data  on 
"Buffalo"  and  "Twin  Otter"  responses  from  other  sources  are  not 
currently  available. 

Accordingly,  it  is  possible  to  say  at  this  time  only  that 
this  report  is  consistent  with  Reference  1  and  that  both  sets  of 
results  are  "reasonable".  The  time  constants,  frequencies,  and 
damping  ratios  of  the  various  modes  presented  in  Appendix  D  of 
Reference  1  agree  with  results  presented  in  this  report.  The 
values  of  these  parameters  are  in  the  expected  ranges ,  and  show 
the  normal  variation  with  airspeed  for  each  aircraft.  Similarly, 
control  power  values  appear  to  be  within  the  expected  ranges  and  in 
proper  proportions. 

Responses  shown  in  this  report  are  for  the  Cruise  Flight  Con¬ 
dition.  For  the  "Buffalo"  this  is  level  flight  at  400  fps  and 


10,000  ft  altitude  with  a  gross  weight  of  40,000  lbs.  For  the 
"Twin  Ottsr",  cruise  is  defined  as  level  flight  at  278  fps  and 
10,000  feet  with  a  gross  weight  of  12,000  lbs. 

Figure  5  shows  the  response  in  pitch  rate  q,  pitch  angle  0, 
angle  of  attack  a,  altitude  rate  li,  and  forward  speed  U  resulting 
from  a  1®  step  elevator  input  6g  for  the  "Buffalo'.*.  Lateral 
degrees  of  freedom  were  suppressed  during  this  run.  This  figure 
compares  v;ith  Figure  D1  of  Reference  1. 

Figure  6  shows  the  same  information  for  the  "Twin  Otter".  This 
figure  corresponds  to  Figure  D13  of  Reference  1. 

Figures  7  and  8  present  lateral  responses  for  the  "Buffalo". 
Here,  longitudinal  modes  are  suppressed.  Figure  7  shows  the  re¬ 
sponse  in  sideslip  angle  3f  roll  rate  P,  roll  angle  4>,  yaw  rate  R, 
and  yaw  angle  f  resulting  from  a  1®  step  aileron  input' 6-^.  '  Fig¬ 
ure  7  compares  with  Figure  D7  of  Reference  1, 

Figure  8  shows  the  response  in  the  same  parameters  resulting 
from  a  1®  step  rudder  input  6^..  This  figure  corresponds  to  Fig¬ 
ure  D8  of  Reference  1. 

Figures  9  and  10  present  lat  "^ral  responses  for  the  "Twin  Otter" 
for  1®  aileron  and  rudder  inputs,  respectively.  These  figures 
correspond  to  Figures  D19  and  D20  of  Reference  1. 
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TABLE  I  , 

SIMULATION  PROGRAM  LISTING' 


•  i 

] 


. .  BLPCX'DATA  '  "■..I.,— . —  .  — ^  .... 

C#***»REVISED  data  for  BurFAUeTER  MG  6/16/71 

REAC’ix>ir/iz - ■  ' 

C0MM0N/C8NST/WE1QHT*RHBSEA#HATMj  AjBiC«SiEPiARl#CTJiCT2*CMT/ 

■  1  CHAtri  CWDALF,  CWQi  CTDCE#  CYB/C^R‘;CYP>  CnB*  CnT^LR/CnRF  I  NiOELCD/ 

S  CNPFIN/CLSiCLRFlNiCLPiCLDLA^TSTAT/ IXi lY/ 1Z/88M  #COF*NVEHiNTYPE/* 

- 3  "•  A'LTBOL - : . . .  “  •  ' 

data  NVEH/ai/  '  ;  t  ,1 

data  WElGHTiRMeS£A,HATM«Al.rBeL./AOOOO,0#.002378,32500.0,.,085/  ' 
data  AiBiC«S/EPlARl/5.2A96,0Al0.i#9A5,0#.0A35/ 
data  Cf1T.CHALFACMDALF/CM0«"ChDLE/0*0A..78i-6.05/-35.6<2.12/'  ' 

_  _ data  CYB«CYR/CYPaCNH«CNDLR.CLB/««362, .368«».055«.101a«107«>-«125/ 

ryXTX  TWmNTDJPTTMTCLRFTTTiCLP#  CLOLA/-  « 169/  •085*  #038/  •  *53#  *20/ 
data  CTl/CT2#OELCD#TSTAT#CDF/,0037#'6,blE-6#*03#22'»00*#.03?/ 
data  IX,IY#I2/273000*#?l5000*#4A7000,/  ,  ' 

end 

BL9CK  DATA  “  ■  ,  i 

C«***»REVtSFO  DATA  FDR  TWIN  BTTER  MG  6''16/71 

- ■RTAL'lX.lYiTr - 

CD“M8N/ceNST/WE!GHT#RHBSEA/HATM»  A#B#C»S|EPI AR1#CT1#CT2#CMTi 

1  CMalF.C«DALF,CHO,CMOLT#CYB,'CYR,CYP#CnB#CnDLR#CNRFIn#DeLCD#  ■ 

2  CNPFIN,CL0iCLRFIN#CLP#CLDLA,TSTAT#IX.lY>IZ#e8M  #CDF#NVFH#NTYPE# 

3  AlFBOU  !  ■■ 

data  NVEH/20/  _  _  _ 

data  re  I  GHT#RH6i>LA#"HATJ1#  AUF  bOL/T200U* '#' *00?378#  32500*0#  •  *092/ 
data  A#B#C#S#EPlARl/5*2#65.0#6*b#A20.0#*0425/ 
data  C!'T#CMALF#CHDAlF#CHQ#CHDlE/0*0*-.78#-6*15#»2‘»*6#1*73/ 
data  CVn#CYR#r.YP#CNB#CNDLR#CUB/»*A92,.A29#-*0RS#*121#*107#-*lC3/ 
data  CNRFIN#CNPFIN#CLRnW#'Cr.P#CLOl.A/-. J68#.*033#*033#-*53#*38/ 
data  CTl#CT2#DELCD#TSTAT#CDF/*00377‘t,9.07E-6#*03b»5750.#,C39/ 

0-ATA  -TT/T  TnZ72¥3(J15‘.7?eDro  *T>fI  OO'O  *  / 

END  !  I 

C  MAIk-  PR0GRAM  :  I  . 

DHENSIQN  0ERIV<12)#VINT(12) 
real  LlFT#LV#LR#LP#LDLA#NV#NR#NPiNDLR, IX# lY# IZ 
X  C8“M9N/aLIPRN/BLPFRQ<ITTB#IBLIP>PRNFHC3#ITTP#lPRN 

■■  C9^5V7T:8T}ST7WET5HT#'RH8SEA,HATH,  '  A,B*C#S#EPIAR.l#CTl#CT2#CMT#  , 

1  CMalF#  CMDALF, CMq, CNDLE#  CYb#  CYR, CYP#  CNB#  CnDLH#  CNHF I  -g, DELCD# 

2  CNPFIN,CLB#CLRFIN,aP;aOLA,TSTAT/Ix#IY#IZ#8nM  #CDF#NvFH#NTYPE# 

3  AlFBOL  i  , 

Ce“'lBN  XSS/YSS#ZSS#RNAX#RNBX#Rnay»PNBY#WNAZ#RN9Z.v,xX#WYY#. 

1  WZZ#XSD#YSD#ZSD#XTAU#YTAU#ZTAU# JX#JY#J2#SGUST#UW#VW#WW 
- C8?t«wr" - THTnR7TnOCE70?TT5r;‘Atrr#vRI#SIGI  .  ‘ 

COmmBN  RHfl#SIG#OYN#CL#CD#YV#YR#YP#LV#CLR,LH#LP#LDLA#NV#CNK/ 

1  NR,CNP/NP#NDLR#CM#0EU#T 

COm-BN  ELTRH#DLE#0LA,DLR#THReT,GDUE#6DUA<GULR,GELTRM# 

1  TRNPNT/SlPBl/PTCHBR  •  I 

C8‘1M0N  ALF#VR>LlFT#ORAG#DALP/THRu5T#UDBT#>o9T#QD8T/VonT#PDOT# 

1  ■■R(5F17T'^^^o■0T;p§ID5T#'PMroBT#XDe^#YD^T#HO0T#U#W#Q#V#P#R,THFTA#PSI# 

2  PHI#X#Y#H#BETA 

COHmBSj  CLO#ALF0O  i  ,  ' 

equivalence  (oERIV(1)#UD8T)# (DERIV(2),wDoT)/ (OERlV(a)#Q[;oT)# 

1  (DERIV(4)#VDBT1/ (DERtV(5)#PDBT)# (OErIV(6)#RD0T1#(DERIV«71#THETDOT 
2)#  (OERIV(8)#PsIDOT)#.(OERIV(9)#PHID0T),  (DERIV«'lO)#XDnT)#  (OER'IVdD# 
-  ■  ^YDOTTnDERTVTIETTHOTTirtVTNTiniUJ/tVINTIBI/WI.IVINTOl/Ql/IVINTIir 
A)#V)# (VINT(5)#P)i (VINT(6)#R)#(VINT(7)/THETA)#{VINT(8)#PGI )# (VINT(9 
5T#Plin#  CVINT{10)#X)#IVINT(in#YI#  (VlNTn2)#H) 
namflist  h#vr,x#y»psi»phi#p#q,r,v#del#ntype#weight 
X  namelist  ELTRM#DLE#OLA#OUR#THR0T  ■  t 

NAmfLIST  TMTHR/TMDLE 
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XXXXXC/'O  XX  ooo 


I 


I 


TABLE  I  (Cont) 


,  namelist  blpfro,prnfrq  .  . 

- namelist  aDLiymiLAiUtllR* IRNPWTTSU^t-VPTCHBRyat.LTRW^ 

NAMel I  ST  XSS*  YSS, ZSS#  XSD#  YSD*  ZSO, XT Au»  YTAU*  ZTAU#  JX/ JY#  JZ 


call  SETP<rr(32,3000«TOOO>OOOOi?nOO) 
call*  SETP9T ( 132/ 3000#  2000i3000#  2000  > 

12  nTYPE'O 
C»f««*  CALL  standby 
S  •  .JMj  031010 

C»*»»t.  3  OEO/SEC/pe INTER  WIDTH  PCS  RATE  OP  TURN 

;RNh  >19. 1 

trtnrr-iyj  ■tJEtSTexctr-WTTJTR' nJR‘  SUTP‘ 

SLPnL«6.73 

c«»»*,*  k  OEQ/BAR-WIdTH  for  pitch 

PTCHBR«1A.325 
C«««*»  2  nEQ/V9LT 

3ELTHM«GDLE«.03A9066 
C»*»*«  2  '3  DES/V9LT  '  '  ’ 

GOLK*. 0116356 

C*«*#»  4/3  OEG/VBLT  * 

GOLAs-. 0232712 

SET  initial  CeNDITIONS#  IDLE  L9PP 

i  H«6C00. 

VR«2C0. 

Y«PSI«PhI»P»O.H«  TMTHR«TMDLE»Z5S*0» 
BLPPROklO. 

PRNPRr.*!* 
jXtjYTTOZ'l  ; 

XSS»14. 

YSS.tlA. 

XS0»2.3 

YS0«1.6  ,  ' 

ZS0«1.  ' 

XTAU«TTAU»ZTAU«t.^ 
iD^La.OS 
call  IP.INITIA 
iNpyTdOS) 
t*tt«  CALL  COMPUTE 
E0M  C31013 
■TEMPTm.TrRa/DEL 
lTTn.TE“P 
TEmp»PHnFR0/0FL 

iTTPiTE^P  T 

IBLIP*IPRN»*1 
80M.32«?/WE1GhT 
00  1  1*1/9 
1  0EHIV(I1«0. 

SIG*EXP(-H/HATM) 

>  RHB.SI.G»RH0SEA  ■' 

CLO»2.*WEIQHT/(RH0»S*VR»VR) 
ALF,30«CL0/A+ALFB0L 
■  “TMrnrrALnior 

IF(sENSESWITCH5)20#21 
•  20  CONTINUE 
U'VR+XSS 
V'YSS  . 

' W'ZSS 
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Q8  TB  32 

— CONTINUE - 

w*v«o. 

’U'VR 

22  CONTINUE 
VWii*Jv««0i 
UW»\/R 

- CD^epr*cL*CL«E.'^T - 

DYN«t5«RHB«VR*vft 

ORAR-CO-OYN'S 

ORARI'ORAG 

AUFI«ALr 

VR!«VR 

■  •“STGT^STS . . 

T«-DEl 

RNA«1..0EL‘/5(tau‘ 

RNAY«lt-DEl/YTAU 
RNAZ.l.-OEL/ZTAU  " 

RNax«XS0«SQRT(2.*0EU/XTAU) 

■  ■'•TITO77VSD»SaRTl2;jn5Ei:7VTAUT — 

HNBZ»ZSD*SgRT(2»«DEU/ZTAU) 

SGUST  'S.-aRT  (  XSS«XSS+YSS»YSS  ) 
call  ARM{0) 

call  enint 

CBNNECT(4O>AER0> 

. . idle  COTP  +  TEST-CARE  CACHE- BF  AeRTT 

10  continue 

11  call  aerb 

idle  LSBp  ♦  FAKE  INT.OO  USINf,  T/F/T.L.aS.L*  IF  INT.  SySTEK  DOaN 
SKS  03000A 
BRU  S«1 

EOT*  03000r  -  -  -  — 

EOn  030200 

IP(SENSESWrTCHl)12#10 

end 

subroutine  aerb 

0lf*CNSl3N  BTE<3»3) 

Ol*TNSlgN-OER-IV{12rA\nNTTI2)  ' 
real  LIFTaLV#LR«LP»UOI.A,nV#nRaNP/NDLR/U»IYaI2 
X  C8M*T0N/?LlPRN/BLPFRn» ITTD/IBLIP/PRnER3/ITTPaIPRN 

C0‘1^‘8N/CBNST/WEIGHT»HMBSEA#HATM,  A/B/CaS»EPIAK1/CT1*CTP/CMT» 

1  CMALFi CNOALF#  CT*a«  CNOLEa  CYB/ CYR*  CYPa  CNDa  CNDLR#  CNRF 1 N/ DELCUi 

2  CNPF1N»CLB#CLHFIN»CLP*CLDLA/TSTAT/ Ix# IY*IZ»BBM  /CDF#Nveh»NTYPL» 

■3-  ‘ALTHOL -  •  - - - 

C6«y8N  XSS/YSS#ZSSiRNAX/RNOX»RNAY*RNBY»RNAZ*RNBZ»wXX»WYY< 

1  N7ZiX3O/YSD/ZS0»XTAJ»YTAU/ZTAl)/JX/ JY»JZ#SGUST»UW*VW#WW 
COt-MBN  THTHR,T‘'DLF.ADRABl,ALfI*VRI/SlGJ 

ce‘'‘*g\'  RH8,SIG»GYN,CL,C0#YV#YR,YP,LV»CLR,LR»LP,LDLA»NViCNR, 

1  NR,cnP»NP,NOlR/C»#OEL»T 

C9*'*3N  ■  FLTRM,OLE/OLA,DLR»TnRBT,3DLE#GDLA,r.DLR,GELTRM, 

1  TRNPNT»SLP0L»PTCHOH 

CB-'-BN  ALF,VR, lift, drag, 0ALF/THRjST#UD6T«W08T»30BT,VD9TiPD0T, 

1  RD?T,TiiETD8T, PSID8T, OHI DOT, XD8T,YDBt, HOOT, U/W,i;,V,P,R,  THLTA,PSI, 

2  PHI, X,Y,H, BETA 
C0H“0N  CL0,ALFBC 

EOin VALENCE- (DERI V(1),UDBT)-,  {0ERIV(2),ViDBT),  (0ERIVt3)»QD0T), 

1  (OERIV(4),VOflT), (DEHIV(5I,PD0T), «DERlV(6)*RO0T)i <0ERIV(7), THE  TOOT 
2)/ tDERIV(8),PSl08T), (DERIVO),PMlD3T), (0ERIV(10>/XDnT), (OERIV(ll), 
3YD0T>,(OERlV(i2),HO8T),tVlNT(l),ul»(VlNT(2),li),{VlNT(3),O),(VlNT(4 
A)/V), (VINT{5),P>,(VINT(fi),R), (vlNT(7),THErA),(VlNT(8)«PSI),  (VINTO 
5) /PHI )/  !VINT(io),X),  (ViNTIli  ),Y),  (VHT(12>,H) 
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i 


c 

CJ»nrT-TT»lTTO"Sn5WXL~;~Srr'  T/y  - - 

s  e9i  030000 

C  RECTANGULAR  INTEGRATION 

T"T*0EL 
00  10  I'lnS 

10  VINT(I)»VINT(I)+OERIV(I)«OEL 

C  •  ttstj-rerf . . 

Ole  erom  -lo  v.  down  to  -hb  v.  up 
. .  Ola  prom  »i6  v,  right  to  *15  v,  left 

C*  >t«  DlR  prom  -30  V.  RIGHT  TO  >.30  V,  LEFT 
ELTRM  prom  -15  V.  DOWN  TO  +15  V*  Up 
ThRBT  prom  -3»2  V»  idle  TO  0  V,  PUlL 
■  CALC  ■AOLTSDiECTRH'.ULEJtJU/DLR/THROT) 
eltrm.qeltrm.fltrh 
ola.gola+ola 
DLR.aOLn*DLR 
0LE«GDLE*0LE+EL,TRM  • 

THRe)T«l.  +  .3l25»THR0T 
C  TOTAL  velocity 

VRSC»UWtUW+VW«VW+WW»WH 

VR«SQRT(VRS0) 

C  FALCjLATE  rOEFPIClFHTS 

G15«CXP(-H/HATM1 
RWe.RHOiJEA+Slo 
lP(T.r,E.THTIIR)Q8  TO  81 
IP  (\VF,H.EQ«?)5IG1.1 

THRf>T«0RAGI»(i,+cTl«VRI+CT2iVfil*vRI » / ( S t G I »TSTAT  ) 

81  lP(T.r,E.TMOLE)Gn  TO 
DLE«0. 

82  CONTINUE 

C  CONVAIR  is  SUPFRCHARGE0»»»US£  GIG.1  POR  THRUST  COMR 

lP(VVEH,£Qt?)si3»l 

THRUST*THR0T«SIQ+TSTAT/(1.+CT1»VR+CT2*VRS0) 

C  OY^A'TC  pressure 
DYni.S»Rm8«VRsQ 
C 

sPHi»siN'(PHn 
SPSI»SIN(PSI ) 

STH.3  IN' (THETA) 

CPhi.COs(PHI  ) 

CPSl^CO.TtPSI) 

^  CTh*C3S(THETA) 

C«****b80V»T0  earth  TRANS  MATRIX 
STCS»ST>l*CPSI 
SSCP«SPSI''CPHI 
sssf»sPr.i*SPHi 
BTE(l<l)«CTH»CPni 
btei  u2)»STcs#sPHi-;;scp 
OTETI ; 3 ) »STCS»CPHI ♦SSGP 
BTE(2,1)«CTH»SPSI 
0TF.(2,2)»SSSP»STH+CPHI»CP5I 
6TF(2ii)»SSCP*STH-SPHI*CPSI 
0TE{3/n»ST!( 

BTE(3/2)»-SPHI*CTH 
BTE(3>3)«-CPHI*CTH 
C  NEW  .;IN3  model 

C  F*S'  21/22/23  O0WN  FQR  GUST  IN  X/V/Z  RtSP* 

C  S»S.%  5  SET  FOR  STCA-DY  STATE 

C 

S  SRS  033006 
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S  BRU  IIS  _ 

CTOT'  Gust  (XGUs«RNAX/ftKI&y«  Jk) 

JX»2 

S  BRU  las 

11  XGUS«0« 

SI2  SKS  030007 

S  0RU  13S  _ 

. TTOt  GUST  (yjjOB/RNAYTRWBYiJy ) 

jY»g 

S  BRU  lAG 

13  y3US»0« 

SI*  St<S  030010 

S  0RU  15S 

. c^ti;  sDsrrzinjS/RKiAi/RNBZJ  jZ) 

jZ»2 
08  Tft  3 
15  ZGUS«0* 

3  IF{SEnSESWITCh5>*;S 

*  ceNTiNue  _ 

WXT»y55*<‘XQU5*XSS^*Y5DS*YS'Sl  /SQUST 
WYY 'YSS* ( XQUS» YSS^YGUS»XSS ) /SGUST 
»(ZZ.2SS+2GUS 
GO  TH  6 

5  nXX.XGUS  . 

wYY«YGUG 

w7r*zGyB  ■  — 

6  ceyiiNUc 

c****.  OnOY  AXIS  VELOCITILS  INCLUOIMG  WHDs 

uW«J-(^XX»BTE(lil)+wYY«3TE(?,l)-w22*iJTC(3il) ) 
VW.V.(WXX»BTEtl,2)4.wYYi*BTr(2i2>.wZZ*BTE«3,2)  1 
WW«a-(  'IXX«BTE(1,3)*wYY»BTE«?/3>»wZZ*0TE«3»3)  ) 

C  ANGUE  Tf  AmCK,tirT/DR)TO"‘  •" 

ALF«ATA\2(WW#uW) 

cl«a«alf+clo 

CO»COr*a»CL«EPlARl 

3S»0YN**G 

LIFT»CU«QS 

ORAr,.co*ns  ... 

c  sideslip 

BETA.ATAN2(VW,UW) 

C  RvS«RVS/2  RVSD»RVsO/2  RV*»WVS0/*  RV*H2»RVSDB/* 
RVS».5»RH9»VR*S 
YV«rv5»CY3 
RVS?.»RVS»B  ■  ■  ■ 

RVa«.5*RVSB 

YR«PV*»CYR 

YP»RV**CYP 

LV«RVS8«CL3 

CL«.CLRFIN>.25*Cl 

RV«(B2.RV*«B 

tR»RV*02»CLR 

L‘’»RV*D?*CLP 

L0LA«RS.B«CLDLA 

NVoKVS3*CNB 

CNR.CNRFIN-.aB.icn-OELCC) 

NR'iRv'tBa^CMR 

CNp.CNPFIN-.25.Ct*(l.-A*EPlAHl) 

NP»RV*B2»CNP 

ndlr«';s»b*ckolr 

C  • 

C  calculate  AEReOY^AIlC  FORCf S 
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CAUF-C8SULF) 

- g»T:rystnr.m:r? - • 

C«*«*«i.9NGITUDINAL  EQUATI9NS 
C«*»**If  NTYPE"2i-  fly  ONLY  LATERAL 
lF(NTYPEtNE*c;90  TO  31 

. ifOOT’iROoT^oaflro. -  ■■■ 

QO  TO  32 

X  X“FBRCt  -  ■ 

31  uDOT»OOM*(THRUST«.DRaQ«CALF*LIFT*SALP-WE1GHT*STH>-Q*W+R*V 
c  “Z"FBRCE  - 

WDOT«00»*(-LIFT»CALF»DRAQ»SALF+WElQHT»CTH*CPHn+0»U-P»V 
C  "AUPHA'OOT  - 


DALF»(W00T-WW«UD0T/UW)«CALF*CALF/UW 

•■'C - n-'HOHENT - ^ - = . .  — ■  ■ 

CM*CMT*CMALF*ALF+C*»5«(CHDALF#DALF+CmO*0)/VR*CMDLE*DLE 

aoor»cw»t}s»c7iY  ■“  —  ■ 
c 

C*»»**LATrRAL  EtJUATlONS' 
c**«»«iF  NTYPC'ii  Fly  only  longitudinal 
“  -TPTNTTPEiWrmxr'Tg  •  32"'  * 

VD0T«PD8T»RD0T«0. 

GO  TO  33 
C  Y  FORCE 

32  VO0r*00M«(YV»VW+YR»R*YP*P*Wr.fGrtT*SPHl  )-R*U*P*W 


L^MfiMENT 

"W»T^V*V»tUR*Rn:P«PTL©tA»t)L  A ) /tx 


C  N  f'OMENT 

ROot*(NV*VW»NR*R^NP»P+NOLR*DLR)/iZ 


eulfr  angle  transformation 

33  THETOnT«Q*CPHI.R«SPHl 
— PSTOBT*  (0»SPHl+Ri'CPfTT7/CTH 
PHIOOT'P+PSIDOT'STH 


•**•*  X0eT»Y00T*HDnT  IN  earth-fixed  COORDinATER 
OO  35  I»l<3 

35  ■  TltRTVll*9T«ri*TF(r;mUF8Tr(I,2).V*BTE(b3)iW 
IF{NTYPE.NE>2)G0  TO  3A 

hoot.thetdot-o. 

PSIOOT'R 
phidbt»p 

'3'?  cOnT  INUt"~'  ““ 

lF(SENSCSWlTCH2)71/72 
71  09«THRaT*10t 

•*»»»0LIPS  FOR  time  On  STRIP  CHART  RECORDER 
X-  •nOCTP’lRLlPn 
X  IPRn-IPHN+1 
X  IF< IBLlP»EQf ITTBJQO  TO  63 
X  o8»-25* 

X  GO  TO  64 
X  63  IOlIP‘0 

-X - - . -  ■" 

X  64  continue 
D1*H**005 

02»VR»tl 

03»THETA»143.25 

DA'0LE*143.25 
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L>  X  X  XX 


1 


TABLE  I  (Cont) 


C 

Ct 

C 


D5»ph1«57«3 

D7«DLA*n*.6 


'I*’’’''?!?!  ^  reR  J.ONQlTUDINAL  eNLY<  2  PAR  LATERAL 


0NLY. 


-7U- 


lP{NTyPE»l>66j 70,45 
"D’S-WDBT - ^ 


06»0«57.3 

D7«-ALr*143.-25 

GO  TO  66 
65  D1«P*57,3 

D2»BETA»286.5 

^15y»!5LH*57»3 


66 

72 

c 


04»ps1*28.65 

CONTINUE  .  ”  . 

call  dal (20/ 01,02*03,04, 05,06,07,08) 

continue  ■  *  ■  “■ 

TURk-.-TRNPNT»7.5*PSID0T  '• 

“snP7--i5;y5i:pgr*B^fA - 

HOei  IN  FT/SEC,  INDICATOR  iN  Pt/MIN 


•6*67  VOL'TS/lOOO  CUIMB,  48*33  VBlTS/1000  DESCEND 


200 

20J 

202' 

208 

209 

210 
211 


lF(HDB.T.GT.O.)ReFC.  ♦AOlJBtHOeT 

lP(HOeT.LE.O.)ReiFC.  .*P98«HDBT  ! 

VR  IN  FTTSftr; - CBNVm“Te'<N’BT5— 

VXTB.VRuSgptSQRTIRHB/ftHBSEA) 

L  ■  !  I 

IP(VKTS-175*)200, 210/210  , 

IP(VKTS-125.)?oJ/269/209  ""  ' 

IF(v<TS-44. )2o2/208/208 

•ATRSPOiO. . -  —  -  ,  , 

GO  TO  2U 

A'lRSPD".  •22«{vkTS-44*  ) 

oO  T8  211  I  . 

AiRSPO*  17.84,092»(VKTS-ia5.)  i 

GO  re  211  (  '  ,  .  . 

atrttpd^  Bs.n-.oHgB.ivsms-iTB;) 
continue  '  • 

AlT**029H“80. 
oIRGVR«-PSI»15, 916667 
I«0rRGYR*.01 

DlRGYW»DIRSYR-l»100.  '  '  ’ 

‘  ■“R8Li:jf95.T»Pm  - - -  , 

If (aBS( ROLL) *GT* 99*9) ROLL'S IGNI 99 *9/ ROLL) 

PITCH'S. *PTCHBR*(THETA*ALF80) 
rPm.3,4THR0T«27. 

OVA  here  .  ‘  •  ! 

CA^L  dal (28/ 

T"  ■  TURn,SLIP/R8FC,AIR5PD,AlT/0IRGYR, roll, PITCH, RPM) 

•••••PRlNT-OjT  AT  FREQUENCY  PRNFRQ 
IF(  JPRM.NE!.  ITTPIGB  TB  62 
IPRS'.O 

61  w'’JTE(10H/l0l)T/U0eT/U,V,W,VR,H/HDeT,ALF*57.3/THETA«b7,3,PH’i«S7 
■T'RST.'57.T/P»5~,j;(J557i3/R»57.3,THRUST,THRnT/LlFT/DRAG/0LE*i)7.3/ 
2  0LA*57.3/0LR »57.3/0ETA»57.3 

lOl  FBR^IATI  J  T»*,|'6,2/3J(,suo8T'«,Fll*A,l2X/6u«8,Fli,4/12X/*V»«/Fll. 

1  12Y/«W**,F11,4/ 11X/*VR«8/Fii,4/6X/9H'8/F1 1.4/9X/9HD0T«*,fi1 .4, 

2  10X/$ALF«»/Fu.4i8Y/*THETA»4,Ft 1.4/ioX/8PHI'*,F11.4//4XitpSI»$ 

3  '^11*^M2X/*P»«/F1i.4/12X/40.*,Fi1.4,  j2X,*R.4,Fji.4,7X,4THRUSTi 


/ 
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OHM  TWIN  OTTCN 

Aiinounocil  in  Augvut  1964.  tho  Twin  Otter  m 
«  STOt  trcntport  powered  bv  two  Prott  A 
Whitney  (UAC)  PT6A*20  turboprop  engtnei. 
Dceign  work  tiegen  m  January  1004,  Comtrue* 
tioii  of  an  initial  batch  of  Twin  Ottera  waaitarted 
in  November  of  the  earre  year  and  the  fint  of 
theeo  flew  on  May  20,  1965. 

At  the  Iwginning  of  1067,  w  total  of  62  Twin 
Otter*  had  bU>u  delivered  or  were  on  order,  with 
options  on  1 1  more.  They  included  eight  for  the 
Chilean  Air  Force,  two  for  Trana<AuelraIien 
Airlinen,  one  for  the  Canailian  Dep^ment  of 
Lands  and  ForeMta,  four  fur  Aeralpi  of  Italy, 
one  for  Northern  Connolidated  Airline*,  and 
othera  for  Fitgnni  Airlinm  and  Air  Wueonain, 
USA.  Production  was  uchediiled  to  ho  at  the 
rate  of  «ix  a  month  through  1967. 

Under  development  for  delivery  in  1968  ia  a 
veraion  of  the  Twin  Otter  with  more  powerful 
(640  eehp)  Pratt  A  Whitney  l'T6A*27  turlioprop 
engimw,  longer  nose  to  provide  more  baggage 
•ipace,  and  AUW  of  12,500  lb  (5.670  kg).  The 


following  data  refer  to  the  current  pro<!uction 

model. 

Tyrs:  Twimturboprop  STOL  tranaport. 

WtHoa:  Braced  high^wing  monoplane,  with  a 
aingle  *treamline.aection  bracing  atrut  on  each 
aide.  Wing  section  NACA  6A  aeriea  mean  line; 
NACA  0016  (modified)  ihickneaa  dutribution, 
A*pect  ratio  10.  Conatant  chord  of  6  A  6  at 
(h&S  in).  Diliedral  2*.  Incidence  2*  30‘. 
No  aweepback.  Altmetal  aafedife  structure. 
All-metal  ailerons  which  also  droop  for  use  as 
flaps.  Double  slotted  sll-metai  full  apsii 
trailing-ed^  flspe.  No  s|>oilera,  Trim-tabs  in 
ailerons.  Pneumatic-boot  de-icing  equipment 
optional. 

Fusclaok.  Conventional  all-inelal  semi  mono* 
coque  safe-lifs  structure. 

Tail  Uhit:  Cantilever  all-metal  stiu^ture  of 
high  strength  iduminium  alloys.  Fin  integrsl 
with  fuselsge.  Fixed-incidence  tailpisne.  Trim- 


tabs  in  rudder  and  port  elevator,  latter  iiitrr- 
connected  with  flaps.  Pneumatic  de  icing 
boots  on  tailplane  leading  edge  optional. 

Lanoiho  Geart  Non-retractable  tricycle  tx^*. 
with  steerable  note-wheel.  Itiitbcr  thuok 
absorption  on  inain  units.  Otro  pneumstic 
nose-wheel  shock-sbtorlwr,  <«oodyrsr  insiii 
wheel  tyres  eiie  1 1-00  x  12,  pressure  32  lbf<tq  in 
(2  25  k|fcirf)  Goodyear  note-wheel  t>re  ti/t- 
8  90  X  12‘50,  pressure  31  Ib/tq  m  (218  IqifrnfL 
Goodrich  hvuraulic  brakes.  Pro\uioii  fur 
alternstivo  float  and  ski  gear. 


Pimi.n  l'iA\T  Two  579  rNlip  Pratt  A  NV)iitiic\ 
(L'.\('|  I*T6.\*‘20  lurlHinrop  engines,  r*cii 

driving  a  llsrtcell  three-bisde  reienible. pitch 
fully-fcalhering  nietsl  pro^ieller,  diameter  8  ft 
0  in  (2  44  ni)  Fuel  in  two  tank*  (8  edit)  under 
cabin  floor:  total  rapacity  919  Imp  gallons 
(4.178  litres)  Two  refuelling  points  on  |Kirt 
side  of  fuselage.  Oil  capacity  2  Imp  gsUoii* 
(9  litres)  per  engine.  Hlectnc  de  icing  s}tt4in 
for  propeller*  and  air-iiitake*  optional. 


49  HlvIllMd  CiJitdt  DH€>t  Twin  Ottir  twln*tur9«pro  traniptil 


Aeoonireodation  I  Two  teats  side-by-sMe  on  flight 
dsok,  Seats  for  13-18  psasengeri  in  msin 
cabin.  Cabin  divided  by  bulkhead  into  main 
passenger  or  freight  compartment  and  baggage 
or  toilet  compartment.  Door  on  each  side  of 
main  cabin,  at  rear,  llaggage  compartments  in 
noss  and  aft  of  cabin,  each  with  upwanl- 
hinged  door  on  port  side, 

SYtTXHa:  Hydraubo  system,  pressure  1,500  lb/ 
SQ  n  (105  k|/crrf).  for  flaps,  brakes  and  nose- 
wheel  steering.  No  pneumatic  system.  One 
200A  stsrter-generator  on  each  engine. 
ExjtCTSONica  AND  KquirMCNT:  Kadio  and  radar 
to  eustoreer’s  specification.  Phndfiying  instru¬ 
mentation  standard, 

Dimbksioks,  xxtbbhal: 

Wingspan  65  ft  Oin  (19  81  m) 

Length  overall  49  ft  6  in  (16  09  m) 

Height  overall  18  ft  7  in  (6-66  m) 

Tsilplane  span  21  ft  Om  (6  40  m) 

Wheel  track  12  ft  6  m  (3’78  m) 

Wheelbase  14  ft  Vm  (4  50m) 


Passenger  door  (port  sid*). 

Height  4  A  2in  (1-27  in) 

Width  1  ft  6  in  (0  76  m) 


Height  tosill  3ft  lOu  (1-17 m) 

Passenger  door  (starboard  side): 

Height  3  ft  9|  in  (1*15  m) 


W|.lth  2  ft  6  in  (0  76  m) 

Height  tosill  3ft  I0tn(l-)7m) 

Hsggage  compartment  dinir  (noAc): 

Heighttosill  3  ft  10  in  (M7  in) 

llsggsge  compartment  door  (port,  rear); 

Height  4  ft  2  m  (1-27  m) 

Width  4  ft  8  H)  (1-42  in) 

Height  to  id)  3ft  10  in  (I  17  in) 

DiMEMSIOMS,  INTERVALl 
Cabin,  excluding  flight  deck,  galley  and  baggage 
or  toilet  compartment, 
l-ength  18  ft  6  in  (5  64  rn) 

Max  width  5  ft  3  m  (l'60  m) 

,M*x  height  4  ft  1 1  in  (1-60  m) 

Floor  arc*  80  2  sq  ft  (7'45  nf) 

Volume  384  cu  ft  (10  87  iif) 

Baggage  compartment  (nusi*)  volume 

22  cu  ft  (u  62  nf ) 

Bsggage  compartment  (rear)  volume 

52  cu  ft  (I  47  Ilf) 

Areas: 

Wtiigs,  gross  420  sq  ft  (39  02  rrf) 

Ailerons  (total)  33-2  tq  ft  (3  08  nf) 

Trafling-^ge  flaps  (total)  112  2  sq  ft  (10  42  rrf) 
Fin  48  0  sq  ft  (4-46  mM 

Iludder,  including  tab  34  0  sq  ft  (3*16  nf) 


Tailpisne  lOd  sq  ft  (9  23  iif) 

Klevators,  including  tab  35  sq  ft  (3-25  iif) 
Weioiitri 

l)aAico(>ers|ing  weight,  inrluiling  pilot  ( 1 70lb>> 
77  kg),  radio  (100  lb«i45  kg)  and  full  ml 

6.170  lb  (2.800  kg) 

Max  payload  (for  100  mile  **  160  km  range) 

4,430  lb  (2,010  Lc) 
MaxT-D  weight  11,679  lb  (5,252  Vx) 

.Max  landing  weight  1 1,000  lb  14,990  if;) 

pERroRMAHrE  (at  max  T-0  weight). 

.Max  cruising  smwd  at  10,000  ft  (3,050  m) 

184  ninh  (297  kndi) 

Kcon  cruising  speed  at  10,000  ft  (3.050  m) 

1.56  mph  (251  knih) 
l.anding  spec«l  64  5  mph  (104  kmh) 

Bate  of  climb  at  S/L  1,550  ft  (472  m)  min 
Service  ceding  25,500  ft  (7.770  rn) 

Service  ceding,  one  engine  out  8.500  ft  (2,590  in) 


T-0  to  60  ft  (15  m): 
8TOI. 


CAR  Pt  3 


I.4nding  from  50  ft  (15  m): 
STOL 


CAR  Pt  3 


1. 120ft  (341  III) 
1.700  ft  (518  n.) 

1.020  ft  (311  in) 
2.160  ft  (658  111) 


Range  with  max  fuel,  30  mm  reserve 

92U  nnlr*  (1,480  km) 


Source;  Reference  3 


FIGURE  1 


OHc-s  lurrALo 

DifferencM  between  the  U8  «nd  C*n*di*n 
>ertions  *re  u  follow* : 

•V*TA.  U8  model,  with  S.850  eehp  Cknerol 
Kleetno  Td4<OB<10  turboproM.  Overall  lenfth 
77  ft  4  n  (29-87  m).  Deeifnation  may  be  eh4ai*<l 
foUowi&f  tranefer  of  reejponiibUity  for  airo^ 
ia  ihJa  oaU|ory  frcm  US  Amy  to  VSAJT. 

00«t18i  Canadian  Defence  Force  model,  with 
3.088  eehp  General  Flectiio  T84/P2  tarboprora. 
Overall  lennh  79  ft  0  in  (34  08  tn).  Otherwiae 
limilar  to  CV<7A,  with  only  amall  difference*  in 
performance. 

WiNOi:  Cantilever  high^wing  monoplane.  Wing 
eection  SACA  84»A4I7‘8  (mod)  at  root, 
HACA  63tA818  (rood)  at  tip.  Aipect  ratio 
9*78.  Chord  It  ft  94  in  (3*89  m)  at  root, 
6  ft  11  n  (1>19  m)  at  tip.  Dihedral  0*  inboard 
of  nncellf*.  8*  outboard.  Incidence  3*  30'. 
Sweepbaek  at  quarter  chord  1*  40'.  Con* 
ventional  faiheafe  roulUapar  structure  of  high* 
strength  aluminium  alloy*,  Full  span  double* 
•lottM  aluminium  alloy  flap*,  outbow  sections 
functionmg  ae  ailerons,  Alurouiium  alloy 
slot'lip  spoiler*,  forward  of  inboard  flaps,  are 
actuated  by  Jarry  Hydraulics  unit.  Spoilers 
coupled  to  manually.operated  ailerons  for 
lateral  control,  uncoupled  for  sypimetrical 
ground  operation.  Electric*lty**ctuated  trim* 
tab  in  starboard  aileron.  Geared  tab  in  each 
aileron.  Ruddor>aileron  interconnect  tab  on 
port  aileron.  Outer  wing  leading-cd|ee  fitted 
with  electrically-controlled  fluan  pneumatic 
rubber  de  ioer  boots. 

Fcssuosi  Fail-saf*  structure  of  high.strength 
aluminium  alloy.  Cargo  floor  supported  oy 
longitudinal  keel  members. 

Tail  Uxrri  Cantilsver  structure  of  high-strength 
aluminium  alloy,  with  flaed-incidsnce  tailpipe 
mounted  at  tip  of  fin.  Ekvator  aerodynamic* 
ally  and  masa-balaneed.  Pore  and  trailing 
•erially*hiAged  rudders  arspowsred  by  taadsra 
Jacks  opsrated  by  two  independent  hydraubo 
^tem*  msnufaetured  by  Jarry  Hyoraulioa. 
'wun-tab  on  port  elevator,  spnng-tab  on 
starboard  slsvator.  Elseirically*oootroUsd 
flush  nneumatic  rubber  de-ioer  boot  on  laU* 
plane  leeding  adge. 


0NM;.liffil9  t«iiidvrHpr9F  iTOL  MtHHy  UMtfti 


l.ANuiso  Gear:  Retractable  tricycle  tyi'e. 
Hydraulic  retraction,  nose  uiui  aft.  main  uiiili 
forward.  Jarry  Hydraulics  olco-pneumatie 
•hock-alisorWrs,  Goodrich  tnsiii  wheels  amt 
tyres,  sue  37  00  y  18  00-13,  pressure  481b  ui 
(3  16  liNi/).  Goodrich  note  ahecis  and  lyres 
site  8  90  y  12  80,  pressure  38  Ib/sq  in  (3  67 
k|/crrf),  Goodrich  muili-disc  brakes. 


Rower  Rlakt:  Tao  General  Electric  T64  turln)* 
prop  engines  (details  under  entries  for  in* 
dividual  versions,  aboiel,  erxh  driiing  a 
Hamilton  Standard  63K60.I3  three-blade  pro* 
peller,  diameter  14  ft  6in  (4  43  m).  Furl  in  one 
integral  tank  in  each  inner  aing.  capacity  833 
Imp  gallons  (3.43d  litres)  and  rublier  hag  tanka 
in  each  outer  aing,  cap^ily  336  Imp  gallons 
(1,837  litres).  Total  fuel  capacity  1,738  Imp 
gallons  (7,900  litres).  Refuelling  points  al>o\e 
wings  Slid  III  side  of  fuselage  for  pressure 
reluellmg.  Total  oil  capacity  10  Imp  gallons 
(48-8  litres). 


‘Dimihsioks,  trriRNAL: 

Wi^  span  96  ft  0  in  (29*26  m) 

LenAh  overall: 

CV*7A  77  ft  4b  (23*67  m) 

CC.116  79  ft  Ott  (24  08  m) 

Height  overall  28  ft  8  n  (6*73  m) 

TailpUne  apan  32  ft  0  n  (9*76  m) 

Wheel  track  30  ft  6  n  (9*29  m) 

Wheelb^  27  ft  1 1  n  (8  80  m) 

Cabin  doors  (each  aide): 

Height  6  ft  6b  (1*68  m) 

Width  2  ft  9  in  (0*84  m) 

Height  to  sill  3  ft  lOm  (M7  m) 

EroerMnoy  eiits  (each  side,  below  wing 
leading-edge)  I 

Height  3  ft  4  b  (1*02  m) 

Width  2ft  3in  (0*66  m) 

Height  to  sill  approx  8  ft  0  n  (1*82  m) 
Rear  cargo  loading  door  and  ramp: 

Height  20  ft  9  b  (6-33  m) 

Width  7  ft  b  n  (2*33  m) 

Height  to  ramp  hinge  3  ft  10  in  (1*17  m) 
Dimxhiions.  intsrmal: 

Cabin,  excluding  flight  deck: 


Length,  cargo  floor 
Max  width 
Max  height 
Floor  area 
Vflume 


31  ft  8 10(9-88  m) 
8ft  9n  (2-67  m) 
8  ft  10m  (2-08  m) 
243  8  sq  ft  (22-03  nf) 
1,718  cu  ft  (48-66  of) 


Areas: 

Wings,  gross  948  sq  ft  (87*8  rrf) 

Ailerons  (total)  39  sq  ft  (3-02  nf) 

Trading-edge  flaps  (tola),  includmg  ailerons) 

280  sq  ft  (26-01  nf) 
Spodcro  (total)  28-2  sq  ft  (2-34  nf) 

Fm  92sqft(888nf) 

Rudder,  including  tab  60  sq  ft  (8*87  nf ) 

Tailplane  I5I'5  sq  ft  (14-07  rif ) 

Elevators,  including  tab  81*5  sq  ft  (7‘57  nf) 
WBiovra  ARD  Loadiroi: 

Operating  weight  empty,  including  3  crew  at 
200  lb  (91  k|)  each,  plus  trapped  fuel  and  oil 
and  full  cargo  handling  equipment 

23,187  Ib  (10.808  i|) 


Max  payload 
Max  T*0  weight 
Max  zero-fuel  weight 
Max  landing  weight 
Max  wing  loa>ling 
Max  power  loading 


13,843  D)  (6.279  kg) 
4l.000Ib  (18,098  kg) 
37.000tb  (10.763  kg) 
39.0001b  (17,690  k|) 
43*4lb/sqft(212kg/nf) 
7*2lb/Mhp  (3-27  kg/eshp) 


PERTORMAHCX  (CV<7A.  at  max  T*0  weight): 
Max  level  epeod  at  10,000  ft  (3.080  ro) 

271  mph  (438  kmh) 
Max  perml*s:ble  diving  epeed 

334  mph  (837  kmh) 

Max  cruising  speed  at  10,000  ft  (3,080  ro) 

271  mph  (438  kmh) 

Econ  cruising  speed  at  10,000  ft  (3,080  ro) 

208  roph  (338  kmh) 
Stalling  speed,  40*  flap*  at  39,000  Ib  (17,690  kg) 
AUW  75  mph  (120  kmh) 

Stalling  speed,  flap*  up  at  max  AUW 

108  mph  (169  kroh) 
Rate  of  climb  at  S/L  1,890  ft  (878  ro)  rain 
Service  celling  30,000  ft  (9,180  m) 

Service  ceiling,  one  engine  out 

14,300  ft  (4,360  ro) 
T*0  run  on  firm  dry  eod  1,040  ft  (317  ro) 
T*0  to  80  ft  (18  m)  from  firm  dry  sod 

1.540  ft  (470  m) 

Landing  from  80  ft  (15  m)  on  firm  dry  a^ 

1.120ft  (342  m) 

Landing  run  on  firm  dry  sod  610  ft  (166  m) 


Source:  Reference  3 


Figure  2 
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L;  EARTH  LOCAL  VERTICAL  COORDINATE  FRAME 
C:  EARTH-AIRCRAFT  CONTROL  COORDINATE  FRAME 
A:  AIRCRAFT  BODY  COORDINATE  FRAME 
EULER  ANGLES 


't'  =  ROTATION  ABOUT  AXIS 
0  =  ROTATION  ABOUT  Y^  AXIS 
4)  =  ROTATION  ABOUT  X^  AXIS 


✓ 


EARTH  HORIZONTAL 
^  PLANE 


Figure  3:  Reference  Coordinate  Frames 
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(a)  At  equilibruim 


(b)  Displaced  from  equilibrium 


Figure  4 ;  Sketches  showing  Relationship  of  A 
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0  10  20  3i 


time  -  sec 

FIG, 7  Response  to  1®  Step  Aileron  Input  (Buffalo,  Cruise) 
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FIG.  10  Response  to  1®  Step  Rudder  Input  (Twin  Otter,  Cruise) 


44 


Appendix;  Development  of  Expression  for  Thrust 

The  longitudinal  force  equation  of  Section  II  includes  the 
total  thrust  force  T.  Since  directly  applicable  data  on  the  pro¬ 
pulsive  system  installation  of  the  "Buffalo"  and  "Twin  Otter"  are 
not  available,  an  expression  for  T  is  developed  here  for  use  in 
the  simulation.  Although  the  expression  is  adequate  for  the  simu¬ 
lation  documented  in  this  report,  it  must  be  considered  an  approx¬ 
imate  one. 

Thrust  developed  by  a  propeller  is 


where  P  is  the  power  supplied  to  the  propeller,  V  is  the  velocity 
of  the  propeller  with  respect  to  the  air,  and  rip  is  the  propeller 
efficiency.  Power  supplied  to  the  propeller  is  expressed  in  this 
report  as 

P  =  a  Pq  ? 

where  a  is  the  atmospheric  density  ratio,  is  the  rated  power  out¬ 
put  of  the  engine  at  sea  level,  and  is  the  pilot's  throttle  de¬ 
flection,  expressed  as  a  fraction 
of  the  deflection  for  rated  power. 

Propeller  efficiency,  rip/  is 
obtained  from  Figure  3-17  of 
Reference  2  (reppoduced  here) 


Fk.i  ui;  3-17.  PioimIK'I  iliiiiciu  v  —ihlr). 
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as  a  function  of  advance  ratio  J  and  power  coefficient  Cp.  By 

definition , 

j  =  60V 
ND 


and 

„  .5P/1000 

^  a(N/1000) ^ (D/10) 5 

where  V  is  in  ft/sec,  N  is  propeller  speed  in  rpm,  D  is  propeller 
diameter  in  feet,  and  P  is  power  in  horsepower  units. 

For  the  "Buffalo"  (Figure  2)  with  its  two  T64-GE-10  engines, 

N  =  1160  rpm,  D  =  14.7  ft,  and  Pq  =  2850  ESPIP/engine ,  so,  at  sea 
level , 

Cp  =  ,137 
or 

C  1/3  =  .515 
P 

Entering  Figure  3-17  at  J/C  1/^  =  2.0  gives  =  ‘79,  This  value 

P  P 

of  J/Cpl/^  corresponds  to  J  -  1.03  or  V  =  293  fps .  Therefore 

(2850)  (550)  , 

T  =  .79  _  =  4220  Ibs/engine 

293 

or,  for  two  engines,  8440  lbs.  Repeating  this  calculation  for 
other  values  of  J/Cpl/3  produces  the  required  thrust  vs  speed 
relationship. 

This  thrust  -  speed  curve  can  be  represented  by  an  equation  of 
the  form 


T 


T 


rated  power, 
sea  level 


static 


V 


R 


By  «curve -fitting  techniques,  it  can  be  established  that,  for  the 


"Buffalo", 
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I 


i 


T 

C 


Static 


=  22400  lbs 


rji^  =  J 00 370  fps 
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I 


Cm  =  6.51x10“°  fps“2 

•*•2  •  ’  .  '  • 

I 

I 


I 

The  process  is  .repeated  for  the  "Twin  Otter"  (Figure  1).  For, this 

•  i 

aircraft  (with  two  PT6A-20  engines),  N  =  2200  rpi^,  D  =  8.5  ft,  and 

‘  ’  { 

Pq  =  652  EiSHP/engihe ._  The  required  constants  are  established  as: 


T 

C 


static 


=  5750  lbs 


=  .00378  fps 


-1 


I 


,  .  I  '  !  i  1 

C^2  9.07x10“®,  fps“^' 

These  values  are  tabulated  in  Section  II  where, simulation  input 
quantities  arc  listed.  .  i  i  ,  . 
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